The  Astrophysical  Journal  Supplement  Series,  169:289-327,  2007  April 
©  2007.  The  American  Astronomical  Society.  All  rights  reserved.  Printed  in  U.S.A. 


AN  ATLAS  OF  [N  n]  AND  [O  hi]  IMAGES  AND  SPECTRA  OF  PLANETARY  NEBULAE 

Arsen  R.  Hajian,1  Steven  M.  Movit,1'2  Denis  Trofimov,1’3  Bruce  Balick,4  Yervant  Terzian,5 

6  7  8  9  10 

Kevin  H.  Knuth,  '  Domhnull  Granquist-Fraser,  Karen  A.  Huyser,  Andre  Jalobeanu, 

Dawn  McIntosh,6  Anne  E.  Jaskot,1’11  Stacy  Palen,12'13  and  Nino  Panagia14’15 
Received  2004  June  30;  accepted  2006  November  26 


ABSTRACT 

We  present  an  atlas  of  Hubble  Space  Telescope  images  and  ground-based,  long-slit,  narrowband  spectra  centered 
on  the  6584  A  line  of  [N  n]  and  the  5007  A  line  of  [O  m].  The  spectra  were  obtained  for  a  variety  of  slit  positions 
across  each  target  (as  shown  on  the  images)  in  an  effort  to  account  for  nonspherical  nebular  geometries  in  a  robust 
manner.  We  have  extended  the  prolate  ellipsoidal  shell  model  originally  devised  by  Aaquist,  Zhang,  and  Kwok  to 
generate  synthetic  images,  as  well  as  long-slit  spectra.  Using  this  model,  we  have  derived  basic  parameters  for  the 
subsample  of  PNe  that  present  ellipsoidal  appearances  and  regular  kinematic  patterns.  We  find  differences  between 
our  parameters  for  the  target  PNe  as  compared  to  those  of  previous  studies,  which  we  attribute  to  increased  spatial 
resolution  for  our  image  data  and  the  inclusion  of  kinematic  data  in  the  model  fits.  The  data  and  analysis  presented  in 
this  paper  can  be  combined  with  detections  of  nebular  angular  expansion  rates  to  determine  precise  distances  to  the 
PN  targets. 

Subject  headings:  atlases  —  planetary  nebulae:  general 


1.  INTRODUCTION 

Planetary  nebula  (PN)  distances  are  commonly  determined 
via  indirect  methods.  A  physical  quantity  is  assumed  to  be  com¬ 
mon  to  and  constant  for  all  PNe,  yielding  fractional  distance  un¬ 
certainties  of  >150%,  and  distances  suitable  only  for  statistical 
surveys.  Despite  the  surprisingly  precise  ground-based  results 
from  astrometric  telescopes  and  careful  calibration  (Harris  et  al. 
1997),  only  about  a  dozen  PNe  yield  accurate  results. 

Without  accurate  distances  to  the  vast  majority  of  individual 
PNe,  some  of  the  most  important  astrophysical  properties  of  PNe 
remain  obscure  despite  many  decades  of  observations.  Among 
these  are  mass  and  luminosity;  estimators  for  both  are  strongly 
distance  dependent,  and  they  are  important  for  several  critical 
projects  including  stellar  evolution  and  mass  loss,  the  PN  lumi¬ 
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nosity  function,  galactic  dynamics  and  energy  balance,  and  the 
cosmological  distance  scale. 

We  have  started  a  program  of  individual  distance  measure¬ 
ments  to  a  collection  of  PNe  using  the  expansion  parallax  method 
(Masson  1986;  Palen  et  al.  2002).  Our  strategy  (Hajian  et  al. 
1993;  Reed  et  al.  1999)  is  to  measure  the  angular  expansion  rate 
of  PNe  by  comparing  at  least  two  epochs  of  image  data  obtained 
with  a  time  separation  of  several  years  between  epochs.  The  ex¬ 
pansion  velocity  of  the  PN  is  determined  using  nebular  spectros¬ 
copy  and  a  spatiokinematic  model  of  the  nebula.  The  resulting 
distance  is  therefore  just  the  quotient  of  the  kinematic  expansion 
rate  and  the  angular  expansion  rate. 

Our  immediate  goal  is  to  obtain  high-dispersion,  narrowband, 
long-slit  spectra  of  target  PNe  with  the  slit  oriented  generally 
along  the  nebular  major  and/or  minor  axis.  Multiple  slit  positions 
were  obtained  in  order  to  more  precisely  determine  the  nebular 
kinematic  geometry.  We  also  have  extended  the  prolate  ellipsoi¬ 
dal  shell  model  by  Aaquist  &  Kwok  (1996)  and  Zhang  &  Kwok 
(1998)  to  generate  synthetic  long-slit  spectra  in  addition  to  syn¬ 
thetic  images  of  the  target  PNe.  We  use  this  model  to  determine 
geometric  parameters  (notably  the  inclination  of  the  major  axis 
to  the  line  of  sight)  for  the  nebulae  so  that  precise  expansion  dis¬ 
tances  can  be  computed  and  systematic  distance  errors  can  be  sig¬ 
nificantly  reduced.  We  present  our  collection  of  spatiokinematic 
data  and  the  modeling  results  in  this  paper  as  an  archive  for  fu¬ 
ture  PN  studies.  We  intend  to  follow  this  work  with  a  paper  in  the 
Astrophysical  Journal  reporting  the  nebular  expansion  distances 
and  resulting  astrophysical  implications. 

2.  DESCRIPTION  OF  THE  SOURCES 
AND  OBSERVATIONS 

2.1.  Images 

During  Hubble  Space  Telescope  (HST)  Observing  Cycles  7, 
8,  and  9,  we  obtained  numerous  images  of  PNe  with  the  WFPC2 
camera  using  the  F502N  filter  to  isolate  light  from  the  [O  hi]  line 
and  the  F658N  filter  for  the  [N  n]  line.  The  images  were  reduced 
using  the  standard  STSDAS  tasks  in  1RAF.  The  multiple  expo¬ 
sures  for  each  filter/target  were  processed  to  remove  cosmic  rays 
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TABLE  1 

A  Summary  of  the  Imaging  Observations 


[Nn] 

[O  in] 

^int 

Tint 

PN 

(s) 

GO/ Data  Set 

(s) 

GO/Data  Set 

(1) 

(2) 

(3) 

(4) 

(5) 

BD  +30  -3639 . 

300 

8390/u5hcl203r 

600 

8390/u5hcl201r 

300 

8390/u5hcl204r 

300 

8390/u5hcl202r 

IC  2448 . 

350 

8390/u5hcl405r 

180 

8390/u5hcl403r 

350 

8390/u5hcl406r 

180 

8390/u5hcl404r 

350 

8390/u5hcl407r 

IC  3568 . 

120 

8390/u5hcl508r 

600 

8390/u5hcl503r 

120 

8390/u5hcl509r 

600 

8390/u5hcl504r 

120 

8390/u5hcl50ar 

600 

8390/u5hcl505r 

230 

8390/u5hcl50br 

600 

8390/u5hcl506r 

230 

8390/u5hcl50cr 

600 

8390/u5hcl507r 

IC  418 . 

200 

8773/u66b  1004m 

200 

8773/u66bl001m 

200 

8773/u66b  1005m 

200 

8773/u66b  1002m 

300 

8773/u66b  1006m 

200 

8773/u66b  1003m 

IC  4663  . 

700 

8773/u66bl203m 

400 

8773/u66bl201m 

700 

8773/u66b  1204m 

400 

8773/u66b  1202m 

NGC  3132 . 

350 

8390/u5hc0103r 

400 

8390/u5hc0101r 

350 

8390/u5hc0104r 

400 

8390/u5hc0102r 

NGC  5882 . 

300 

8390/u5hc0305r 

860 

8390/u5hc0304r 

300 

8390/u5hc0306r 

300 

8390/u5hc0307r 

300 

8390/u5hc0308r 

NGC  5979 . 

600 

8390/u5hc0405r 

230 

8390/u5hc0403r 

600 

8390/u5hc0406r 

230 

8390/u5hc0404r 

NGC  6565 . 

400 

8773/u66b0104r 

260 

8773/u66b0101r 

400 

8773/u66b0105m 

260 

8773/u66b0102r 

400 

8773/u66b0106r 

200 

8773/u66b0103r 

NGC  6578 . 

400 

8390/u5hc0605r 

230 

8390/u5hc0603r 

400 

8390/u5hc0606r 

230 

8390/u5hc0604r 

NGC  6741 . 

300 

8773/u66b0404r 

300 

8773/u66b0401r 

260 

8773/u66b0405r 

300 

8773/u66b0402r 

400 

8773/u66b0406r 

200 

8773/u66b0403r 

NGC  6818 . 

500 

8773/u66b0204r 

160 

8773/u66b0201r 

500 

8773/u66b0205r 

160 

8773/u66b0202r 

300 

8773/u66b0206r 

200 

8773/u66b0203r 

NGC  6826 . 

400 

8390/u5hc0702r 

400 

8390/u5hc0703r 

400 

8390/u5hc0704r 

400 

8390/u5hc0705r 

NGC  6884 . 

400 

8390/u5hc0805r 

260 

8390/u5hc0803r 

400 

8390/u5hc0806r 

300 

8390/u5hc0804r 

NGC  6886 . 

200 

8773/u66b0604r 

200 

8773/u66b0601r 

200 

8773/u66b0605r 

200 

8773/u66b0602r 

300 

8773/u66b0606r 

200 

8773/u66b0603r 

NGC  6891 . 

60 

8390/u5hc0907r 

60 

8390/u5hc0902r 

180 

8390/u5hc0908r 

60 

8390/u5hc0904r 

200 

8390/u5hc0909r 

120 

8390/u5hc0905r 

200 

8390/u5hc090ar 

120 

8390/u5hc0906m 

NGC  7026 . 

400 

8773/u66b0704r 

200 

8773/u66b0701r 

400 

8773/u66b0705r 

160 

8773/u66b0702r 

300 

8773/u66b0706r 

160 

8773/u66b0703r 

NGC  7354 . 

400 

8773/u66b0801r 

300 

8773/u66b0802r 

Note. — All  of  the  data  shown  in  this  paper  were  acquired  from  projects  GO 
8390  and  GO  8773  (PI:  Hajian). 


TABLE  2 

A  Summary  of  the  Spectroscopic  Observations 


SLITs/7;nt  (s) 


PN 

(1) 

Obs./Year 

(2) 

[Nn] 

(3) 

[O  III] 

(4) 

BD  +30  -3639 .... 

K99 

2/15 

4/60 

IC  2448 . 

C99 

6/3  x  30 

IC  3568 . 

K99 

1/3  x  240,  1/120 

4/120 

IC  418 . 

C99 

7/20 

7/3  x  20 

IC  4663  . 

C99 

2/120 

NGC  3132 . 

C99 

4/60 

5/3  x  240 

NGC  5882 . 

C99 

12/60 

4/60,  2/3  x  10,  1/10 

NGC  5979 . 

C99 

6/60 

NGC  6565 . 

C99 

8/60 

6/60 

NGC  6578 . 

C99 

9/4  x  30,  1/3  x  120,  4/30 

COO 

3/300 

NGC  6741 . 

C99 

2/30 

6/30 

NGC  6818 . 

C99 

10/120 

NGC  6826 . 

K99 

11/120 

NGC  6884 . 

K99 

2/120 

4/60 

NGC  6886 . 

K99 

4/60 

4/60 

NGC  6891 . 

K99 

6/120 

NGC  7026 . 

K99 

17/various 

12/various 

NGC  7354 . 

K99 

1/120,  5/18 

Note. — K99  =  KPNO  1999,  K00  =  KPNO  2000,  C99  =  CTIO  1999,  and 
COO  =  CTIO  2000. 

using  the  IRAF  task  crrej,  and  the  distortion  correction  was 
implemented  using  drizzle  (Fruchter  &  Flook  2002)  and  the 
Trauger  coefficients  (Holtzman  et  al.  1995).  Leakage  of  the  Ha 
line  into  the  sharp  edges  of  the  F658N  bandpass  is  unavoidable 
due  to  the  4%-5%  transmission  of  the  filter  at  656.3  nm.  Correct¬ 
ing  for  this  effect  is  virtually  impossible  for  this  study,  since  even 
the  Doppler  shifts  of  the  Ha  gas  can  locally  change  the  transmis¬ 
sion  through  the  F658N  filter.  However,  since  the  Ha  and  [N  n] 
morphologies  are  usually  drastically  different  in  a  particular  PN 
(especially  for  high-ionization  PNe),  the  effects  of  Ha  leakage 
do  not  affect  the  analysis  in  this  paper. 

A  summary  of  the  image  observing  parameters  and  sources 
is  shown  in  Table  1.  Column  (1)  denotes  the  nebula,  and  col¬ 
umns  (2)-(5)  show  the  integration  times  and  references  for  the 
[N  ii]  and  [O  hi]  images,  respectively.  For  all  but  one  of  the  PNe, 
the  image  used  corresponds  to  the  1 1ST  WFPC2  image  of  the  PN 
in  the  light  ofthe  emission  line  corresponding  to  the  long-slit  spec¬ 
tra.  The  HST  data  sets  that  were  used  in  constructing  the  images 
shown  in  the  paper  are  tabulated  in  Table  2.  All  of  the  images  in 
this  compilation  were  obtained  by  us  as  part  of  HST  GO  8773 
and  8390  (these  images  were  duplicated  in  GO  7501).  In  all  of 
the  figures,  north  is  toward  the  top  of  the  page  and  east  toward  the 
left-hand  side. 

2.2.  Spectroscopy 

Echelle  data  from  a  target  list  of  PNe  were  obtained  using  the 
4  m  telescopes  at  the  Kitt  Peak  National  Observatory  (KPNO) 
and  the  Cerro  Tololo  Inter- American  Observatory  (CTIO).  The 
KPNO  observations  were  conducted  on  the  nights  of  1999  June 
25-28  and  2000  April  18-19,  and  the  CTIO  observations  were 
conducted  on  the  nights  of  1999  May  28-29  and  2000  June  17 
and  20. 

The  spectroscopic  observations  were  obtained  in  the  light  of 
[N  n]  (6584  A)  and  [O  hi]  (5007  A)  whenever  possible  for  three 
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Fig.  1. — BD  +30  —3639  [N  n]  images  from  HST and  spectra  from  KPNO  1999. 


of  the  four  observing  runs.  For  the  run  at  CTIO  during  2000, 
inclement  weather  permitted  only  the  acquisition  of  [N  n]  data. 
Seeing  was  variable  throughout  the  observations  on  timescales 
of  a  few  hours  and  generally  averaged  between  1 .0"  and  1 .5".  The 
atmospheric  conditions  were  not  cmcial  to  the  observations  due  to 
the  narrow  observing  bandpass,  the  large  aperture,  and  the  gener¬ 
ally  high  surface  brightness  of  the  target  nebulae.  For  all  four 
observing  runs,  the  80"  decker  and  1.5"  slit  were  used,  and  the 
echelle  cross-disperser  was  replaced  with  a  mirror,  so  that  only 
one  order  was  obtained  centered  near  the  wavelengths  of  inter¬ 
est.  The  long-slit  capability  of  the  spectrometer  required  the  use 
of  an  order-separation  filter.  To  maximize  the  observations,  either 
[N  ii]  or  [O  hi]  data  were  obtained  so  as  to  avoid  switching  the  grat¬ 
ing  and  recalibrating  the  instrument  during  the  evening.  The  filters 
had  total  bandwidths  of  «45  A  for  the  [N  n]  observations  and  80  A 
for  the  [O  hi]  observations.  The  instrumental  dispersion  for  the 
[N  n]  observations  was  0.081, 0.083, 0.082,  and  0.08 1  A  pixel' 1 
for  the  KPNO  1999,  KPNO  2000,  CTIO  1999,  and  CTIO  2000 
data  sets,  respectively,  corresponding  to  «3.7  km  s'1.  For  the 


[O  hi]  observations,  the  instrument  dispersion  was  0.062,  0.062, 
and  0.060  A  pixel'1  for  the  KPNO  1999,  KPNO  2000,  and  CTIO 
1999  data  sets,  respectively,  also  corresponding  to  ~3.7  km  s'1. 
By  carefully  examining  unresolved  features,  we  estimate  that  the 
spectral  resolution  of  the  observations  ranges  between  1.8  and 
2.3  pixels.  We  expect  that  systematic  errors  inherent  to  the  spec¬ 
trometers  and  the  wavelength  calibration  process  will  limit  us  to 
sal  km  s'1  accuracy.  In  the  spatial  dimension  we  applied  2  pixel 
on-chip  binning  for  both  CTIO  data  sets  and  the  KPNO  1999 
data  set,  resulting  in  0.50"  and  0.42"  pixel'1,  respectively.  For 
the  KPNO  1999  data  set,  no  binning  was  applied,  and  the  result¬ 
ing  dispersion  is  0.21"  pixel'1.  The  spectroscopic  observations 
are  summarized  in  Table  2.  Column  (2)  shows  the  observing  run 
that  resulted  in  the  data  (C  =  CTIO  and  K  =  KPNO),  and  columns 
(3)  and  (4)  show  the  number  of  slit  positions  and  corresponding 
integration  times  for  the  [N  n]  and  [O  hi]  images,  respectively. 

All  of  the  PNe  studied  in  this  paper  are  spatially  resolved  from 
the  ground.  Each  PN  is  spanned  by  at  least  five  spatial  resolution 
elements,  so  the  finite  width  of  the  slit,  blurring  of  the  apparent 


Fig.  2. — BD  +30  —3639  [N  ii]  model  image  and  spectra. 
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slit  position  on  the  sky  due  to  seeing,  telescope  guiding,  and 
differential  atmospheric  refraction  are  minimal. 

We  bias-subtracted  the  echelle  spectra  and  performed  flat- 
fielding  using  dome  flats.  We  used  thorium-argon  lamp  spectra  to 
solve  for  geometric  distortions  by  iteratively  executing  the  1RAF 
tasks  identify,  reidentify,  fitcoords,  and  transform.  In 
some  instances,  we  used  sky  lines  in  the  images  to  test  the  ac¬ 
curacy  of  these  corrections.  There  were  no  significant  deviations 
at  the  0.1  pixel  level.  We  also  used  the  lamp  spectra  to  calibrate 
the  wavelength  scale.  The  resulting  two-dimensional  spectra  are 
shown  in  Figures  1-56,  along  with  an  image  of  the  PN  with  the 
slit  positions  located.  In  order  to  determine  the  slit  orientation 
(i.e.,  to  resolve  the  180°  ambiguity  in  the  slit  position  angle),  the 
letters  denoting  the  slit  locations  on  the  HST  images  correspond 
to  the  end  of  the  slit  identified  in  the  spectra. 

3.  ANALYSIS 

3.1.  Description  of  Spatiokinematic  Model 

The  goal  of  the  modeling  described  below  is  to  infer  the 
spatiokinematic  structure  of  a  PN  at  a  particular  epoch  by  fitting 
models  to  image  and  long-slit  spectral  data.  Without  this  infor¬ 
mation,  large  errors  can  be  incurred  when  radial  expansion  ve¬ 
locities  are  married  to  plane-of-the-sky  expansions  in  order  to 
infer  distances.  For  the  subsample  of  PNe  described  above,  the 
nebular  shells  are  acceptably  approximated  by  projected,  prolate, 
ellipsoidal  shells.  We  adopt  this  formalism,  which  is  known  as 
the  prolate  ellipsoidal  shell  (PES)  model  (Aaquist  &  Kwok  1996; 
Masson  1989a,  1989b;  Zhang  &  Kwok  1998),  and  sketch  our  im¬ 
plementation  of  the  classical  model  below. 

Following  the  derivation  of  Aaquist  &  Kwok  ( 1 996,  Appendix), 
we  begin  with  a  central  star  that  is  a  source  of  ionizing  photons  and 
is  therefore  surrounded  by  an  ionized  hydrogen  density  distribu¬ 
tion,  n(9,  </>,  r ).  We  treat  this  situation  similar  to  that  of  a  Stromgren 
sphere:  the  star  ionizes  the  gas  to  a  radius  Rq(9,  f)  depending  on 
the  intervening  column  density  of  hydrogen  gas  (i.e.,  the  nebula 
is  ionization  bounded).  In  addition,  we  assume  that  the  inner  cav¬ 
ity  of  the  PN  out  to  a  distance  Rf9,  <j>)  has  been  vacated  by  the 
high-velocity  wind  from  the  central  star.  In  addition,  we  have 
assumed  that  all  targets  are  optically  thin  at  the  observed  wave¬ 
lengths,  which  may  not  be  correct  in  some  cases.  If  we  equate  the 


number  of  ionizing  photons  along  a  line  of  sight  ( left-hand  side 
below)  to  the  number  of  local  recombinations  (right-hand  side 
below),  we  arrive  at 


dL 

dLl 


agn  (9,  </>,  r)r2  dr, 


(1) 


where  Q  is  the  solid  angle,  L  is  the  stellar  ionizing  luminosity,  and 
as  is  the  hydrogen  recombination  coefficient  (excluding  captures 
to  the  ground  state). 

Now  we  separate  the  density  distribution  into  three  components: 

n(9,  r)  =  (2) 

where  n0  =  «(0, 0,  Rf)  [i.e.,  r)r(r  =  Rt)  =  I ,  rjo(9  =  0)  =  1, 
=  0)  =  1],  the  angle  6  is  with  respect  to  the  polar  axis,  and 
the  angle  <j>  is  the  azimuthal  coordinate.  Note  that  the  model  PN 
is  symmetric  with  respect  to  the  angle  f.  If  we  choose  a  radial 
power-law  distribution  of  the  form 


then  equation  (1)  has  an  analytic  solution  of  the  form 

dL  _  aBn$ri$riffi 
dfl  ~  -2  7  +  3 


n  \-2l+3 

-1 

Rt 


(4) 


The  above  model  is  physically  meaningful  only  if  dL/(dfl)  >  0 
and  if  R ,  <  Ro,  which  can  only  be  satisfied  for  7  <  1.5.  The 
azimuthal  density  distribution  is  parameterized  by 


99  = 


(1 

(1 


P) 

P) 


(2 9\a  (3  ' 

w 

(2n  -  2 9\a  (3  ' 

v- ^ 


for  0  <  9  <  ^ , 

-  ~  2 

for  —  <  9  <  tt, 

2  ~ 

(5) 


where  a  is  the  latitudinal  density  gradient  and  / 3  is  the  ratio  of  the 
polar  to  equatorial  density. 


Fig.  5. —  IC  2448  [O  m]  images  from  HST and  spectra  from  CTIO  1999. 


Fig.  6 a 


Fig.  6. — (a)  IC  2448  [0  iii]  core  model  image  and  spectra  of  nebular  core.  ( b )  IC  2448  [0  iii]  halo  model  image  and  spectra  of  nebular  halo. 


294 


Fig.  6b 
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Fig.  7. — IC  3568  [N  n]  images  from  HST and  spectra  from  KPNO  1999. 


Fig.  8. — IC  3568  [N  n]  model  image  and  spectra. 


Fig.  9. — IC  3568  [O  iii]  images  from/KT and  spectra  from  KPNO  1999. 
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Fig.  10  b 


Fig.  10. —  (a)  IC  3568  [O  hi]  core  model  image  and  spectra  of  the  nebular  core.  ( b )  IC  3568  [O  hi]  halo  model  image  and  spectra  of  the  nebular  halo. 
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Fig.  11. —  IC  418  [N  ii]  images  from  HST and  spectra  from  CTIO  1999. 
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Fig.  15. — IC  4663  [N  n]  images  from  HST and  spectra  from  CTIO  1999. 
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Fig.  21. — NGC  5882  [N  n]  images  from  HST and  spectra  from  CTIO  1999. 


A 

B 

c 

D 

o 

V/ 

G 

, ,  100  km/s 

E 

F 

H  , _ „ 

• 

X 

X 

X 

p 

1 

J 

K 

L 

X 

o 

o 

Fig.  22.— NGC 

5882  [N  ill  m 

odel  image  and  spectra. 

s 


B 


C 


Fig.  23. — NGC  5882  [O  m]  images  from  HST and  spectra  from  CTIO  1999. 
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Fig.  25. — NGC  5979  [0  iii]  images  from  HST and  spectra  from  CTIO  1999. 


Fig.  27. — NGC  6565  [N  n]  images  from  HST and  spectra  from  CTIO  1999. 


Fig.  28. — NGC  6565  [N  n]  model  image  and  spectra. 


Fig.  29. — NGC  6565  [O  hi]  images  from  HST and  spectra  from  CTIO  1999. 


Fig.  30. — NGC  6565  [O  hi]  model  image  and  spectra. 
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Fig.  33. — NGC  6578  [O  hi]  images  from  HST and  spectra  from  CTIO  1999. 


Fig.  34. — NGC  6578  [O  hi]  model  image  and  spectra. 
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Fig.  41. — NGC  6826  [N  ii]  images  from  HST and  spectra  from  KPNO  1999. 


For  the  kinematic  portion  of  the  model,  we  assumed  that  since 
all  of  the  nebular  gas  was  present  at  the  location  of  the  central  star 
at  some  point  in  time,  the  expansion  velocity  of  the  gas  is  purely 
radial  and  the  magnitude  of  the  velocity  is  proportional  to  the 
distance  from  the  central  star. 

3.2.  Implementation  of  Fitting  and  Parameter  Estimation 

To  implement  the  model  described  above,  we  begin  with  an 
ellipsoidal  surface,  /?,  ,  in  the  frame  of  the  PN  defined  by  two  pa¬ 
rameters:  the  length  of  the  major  axis,  R,t  maj ,  and  the  axial  ratio, 
Ar.  The  third  parameter  is  the  ionization  parameter,  f: 

T  _  dL/{dn) 

^°  /  \2  ’  W 

(«0 vt)  aB 

which  was  used  in  conjunction  with  a,  (3,  and  7  to  describe  the 
density  distribution  and  then  to  compute  Rq(0 ,  </>).  We  set  the  in¬ 
tensity  at  each  voxel  (volume  element)  equal  to  aRn2  provided 
the  voxel  is  between  R,(0,  <j>)  and  Rq(9,  </>). 

We  then  rotate  the  three-dimensional  structure  of  the  PN  into 
the  observer’s  frame  using  the  angles  i  and  O',  which  are  the  in¬ 
clination  relative  to  the  line  of  sight  and  the  position  angle  on  the 
sky.  The  intensities  are  integrated  along  the  line  of  sight  to  form 


a  two-dimensional  model  representing  the  brightness  of  the  PN 
on  the  sky.  Kinematic  modeling  requires  the  addition  of  a  single 
parameter:  the  equatorial  velocity,  vq  =  v(n/2, 0 ,/?,■).  The  veloc¬ 
ity  in  each  voxel  is  oriented  away  from  the  central  star,  and  the 
magnitude  is  proportional  to  the  distance  from  the  central  star 
(normalized  to  i>o).  After  rotating  the  PN  to  the  desired  orienta¬ 
tion,  model  long-slit  spectra  can  be  synthesized  for  each  image 
pixel  by  integrating  the  intensity-weighted  radial  velocity  for  each 
voxel  in  the  line  of  sight.  In  summary,  an  image  model  is  uniquely 
defined  by  a  selection  of/o,  a,  /?,  7,  AR,  i,  and  O',  and  the 
kinematic  model  requires  vo  in  addition. 

We  generated  synthetic  images  and  spectra  at  the  slit  positions 
corresponding  to  the  observed  nebular  major  and  minor  axes. 
Since  most  of  our  PNe  have  thin  shells,  we  adopted  7  =  0  for 
all  model  realizations  (see  Zhang  &  Kwok  1998).  To  derive  the 
“best-fit”  model  parameters,  matches  between  model  and  data 
images/spectra  were  performed  by  eye.  This  method  works  well 
due  to  the  orthogonality  of  the  spatial  and  spectral  data  sets;  the 
images  constrain  primarily  the  nebular  structure  orthogonal  to  the 
line  of  sight,  while  the  spectra  constrain  nebular  structures  primar¬ 
ily  along  the  line  of  sight.  Although  far  from  an  optimal  technique 
for  matching  models  to  data,  Zhang  &  Kwok  (1998)  did  the  fits 
by  eye  as  well.  It  is  unclear  whether  a  numerical  technique  for 
searching  parameter  space  for  the  model  that  best  matches  the 
observed  data  is  useful,  since  the  PES  model  produces  the  most 
reasonable  results  in  cases  of  approximately  ellipsoidal  shells 
(e.g.,  the  relationship  between  radial  velocity  and  radial  structure 
is  constrained  by  the  ellipsoidal  geometry). 

For  PNe  with  model  fits,  we  display  the  best-fit  model  images 
and  spectra  following  the  data  in  Figures  1-56.  Seven  of  the  nine 
best-fit  model  parameters  are  listed  in  Table  3:  we  do  not  list  R , 
since  the  scale  is  arbitrary,  and  7  was  not  permitted  to  vary.  Best- 
fit  models  were  not  made  for  some  of  the  observed  PNe  due  to 
the  presence  of  complicated  slit  orientations,  unresolved  kine¬ 
matic  structure,  and  irregular  nebular  structure.  Uncertainties  in 
the  fitted  parameters  were  crudely  estimated  based  on  the  quality 
of  the  model  fits.  To  be  specific,  uncertainties  were  chosen  by 
noting  the  positions  in  parameter  space  where  the  solution  was 
an  obvious  mismatch  to  the  data.  In  this  way,  we  are  probably 
overestimating  the  true  errors  in  the  parameters  and  have  erred 
on  the  side  of  caution.  We  have  assigned  uncertainties  to  the  axial 
ratios,  inclination  angles,  and  expansion  velocities  of  0.15,  10°, 
and  2  km  s'1,  respectively.  For  those  targets  analyzed  by  Zhang 
&  Kwok  (1998)  and  in  this  work,  we  tabulate  their  best-fit  pa¬ 
rameters.  The  results  from  Zhang  &  Kwok  (1998)  were  computed 
for  the  Ha  line,  and  as  such  are  generally  applicable  to  the  [O  hi] 
data  sets  in  this  paper.  The  exceptions  are  low-ionization,  young 
PNe  such  as  IC  418  or  BD  +30  —3639.  The  agreement  between 
the  Ha  and  [N  11]  data  is  expected  to  be  somewhat  worse,  due 
mainly  to  the  difference  in  ionization  state  between  the  two  ions. 
Since  [N  11]  has  a  significantly  lower  ionization  potential  than 
[O  hi],  we  expect  the  [N 11]  emission  to  be  farther  away  from  the 
central  stars,  and  this  effect  is  visible  in  our  results  in  most 
cases.  In  all  of  the  cases  examined  in  this  paper,  we  find  that  the 
[N  11]  gas  has  a  larger  apparent  expansion  velocity  than  the  [0 111] 
gas.  This  inference  is  expected  due  to  the  lower  ionization  poten¬ 
tial  of  the  [  N  11]  and  the  generally  larger  distance  of  this  gas  from 
the  central  star. 

The  main  differences  between  the  Zhang  &  Kwok  (1998) 
study  and  our  analysis  are  due  to  ( 1 )  our  use  of  higher  spatial  res¬ 
olution  image  data  (i.e.,  WFPC2  vs.  ground-based)  and  (2)  our 
inclusion  of  kinematic  information  in  the  PES  model.  The  im¬ 
proved  spatial  resolution  of  our  data  results  in  more  precise  (and 
different)  estimates  of  the  shell  thickness  and  in  some  cases,  the 


Fig.  42 a  Fig.  42 b 


Fig.  42. — (a)  NGC  6826  [N  n]  core  model  image  and  spectra  of  the  nebular  core.  ( b )  NGC  6826  [N  n]  halo  model  image  and  spectra  of  the  nebular  halo. 


Fig.  44. — NGC  6884  [O  hi]  model  images  and  spectra. 
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Fig.  45. — NGC  6886  [N  n]  images  from  HST and  spectra  from  KPNO  1999. 


Fig.  46. — NGC  6886  [N  n]  model  image  and  spectra. 
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Fig.  49. — NGC  6891  [O  m]  images  from  1 1ST  and  spectra  from  KPNO  1999. 


A 


B 


C 


D 


Fig.  50. — NGC  6891  [0  hi]  model  image  and  spectra. 


Fig.  51. — NGC  7026  [N  n]  images  from  HST and  spectra  from  KPNO  1999. 
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Fig.  53. — NGC  7026  [O  in]  images  from  HST and  spectra  from  KPNO  1999. 
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Fig.  54. — NGC  7026  [O  hi]  model  image  and  spectra. 
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Fig.  55. — NGC  7354  [O  in]  images  from  HST and  spectra  from  KPNO  1999. 
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Fig.  56. — NGC  7354  [O  hi]  model  image  and  spectra. 
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TABLE  3 

Tabular  Index  of  Model  Fits  to  Spectroscopy  and  Imaging  Observations 


PN 

Species 

Data  Figure 

Model  Figure 

a 

P 

Ar 

h 

i 

(deg) 

9' 

(deg) 

VQ 

(km  s-1) 

BD  +30  -3639 . 

[N  ii] 

i 

2 

0.40 

0.15 

0.45 

12.50 

28.6 

96.7 

11.0 

[O  m] 

3 

4 

0.40 

0.15 

0.40 

12.50 

28.6 

96.7 

9.0 

Ho 

ZK 

ZK 

1.00 

0.90 

0.60 

30.0 

IC  2448  (core) . 

[O  m] 

5 

6  a 

0.50 

0.85 

0.85 

5.00 

85.9 

135.0 

12.0 

IC  2448  (halo) . 

[O  m] 

5 

6  * 

0.50 

0.85 

0.84 

4.00 

85.9 

135.0 

22.0 

IC  2448 . 

Ho 

ZK 

ZK 

1.00 

0.90 

0.67 

30.0 

IC  3568 . 

[N  ii] 

7 

8 

1.00 

0.85 

1.00 

3.00 

90.0 

0.0 

14.0 

IC  3568  (core) . 

[O  m] 

9 

10a 

1.00 

0.85 

1.00 

3.00 

90.0 

0.0 

12.0 

IC  3568  (halo) . 

[O  m] 

9 

10* 

1.00 

0.00 

0.98 

3.00 

90.0 

0.0 

21.4 

IC  3568 . 

Ho 

ZK 

ZK 

1.00 

0.50 

1.00 

30.0 

IC  418 . 

[N  ii] 

11 

12 

0.30 

0.30 

0.95 

10.00 

57.3 

156.9 

13.0 

[O  m] 

13 

14 

0.30 

0.30 

0.82 

7.50 

57.3 

156.9 

10.0 

Ho 

ZK 

ZK 

1.00 

0.50 

0.83 

30.0 

IC  4663  . 

[N  ii] 

15 

16 

2.00 

0.50 

0.62 

5.00 

97.4 

83.2 

30.0 

Ho 

ZK 

ZK 

1.00 

0.30 

0.75 

32.0 

NGC  3132 . 

[N  ii] 

17 

18 

0.60 

0.50 

0.58 

5.00 

122.6 

150.0 

33.0 

[O  m] 

19 

20 

0.60 

0.50 

0.58 

7.00 

122.6 

150.0 

14.0 

Ho 

ZK 

ZK 

1.00 

0.80 

0.71 

70.0 

NGC  5882 . 

[N  ii] 

21 

22 

0.75 

0.50 

0.58 

4.00 

75.7 

0.0 

26.0 

[O  m] 

23 

24 

0.75 

0.50 

0.58 

3.50 

75.7 

0.0 

24.0 

NGC  5979 . 

[O  m] 

25 

26 

0.75 

0.70 

0.65 

2.50 

57.3 

152.2 

18.0 

NGC  6565 . 

[N  ii] 

27 

28 

1.00 

0.40 

0.36 

3.25 

166.2 

0.0 

17.5 

[O  m] 

29 

30 

1.00 

0.40 

0.36 

4.75 

166.2 

0.0 

15.0 

Ho 

ZK 

ZK 

1.00 

0.90 

0.67 

40.0 

NGC  6578 . 

[N  ii] 

31 

32 

1.00 

0.50 

0.56 

2.75 

157.0 

123.0 

16.0 

[O  m] 

33 

34 

1.00 

0.50 

0.56 

2.75 

157.0 

133.0 

15.0 

NGC  6741 . 

[N  ii] 

35 

36 

1.50 

0.50 

0.54 

5.00 

120.3 

104.0 

26.0 

[O  m] 

37 

38 

1.50 

0.50 

0.54 

6.00 

120.3 

104.0 

20.0 

NGC  6818 . 

[N  ii] 

39 

40 

1.50 

0.70 

0.76 

2.50 

65.9 

13.0 

30.0 

NGC  6826  (core) . 

[N  ii] 

41 

42a 

0.75 

0.50 

0.72 

5.00 

118.0 

120.0 

10.0 

NGC  6826  (halo) . 

[N  ii] 

41 

42* 

0.75 

0.50 

0.79 

3.50 

118.0 

120.0 

22.0 

NGC  6884 . 

[O  m] 

43 

44 

0.50 

0.50 

0.62 

4.50 

117.5 

20.0 

15.0 

NGC  6886 . 

[N  ii] 

45 

46 

0.30 

0.50 

0.34 

12.00 

90.0 

130.0 

17.5 

[O  m] 

47 

48 

0.30 

0.50 

0.34 

12.00 

90.0 

130.0 

13.0 

NGC  6891 . 

[O  m] 

49 

50 

1.00 

0.35 

0.66 

3.00 

120.3 

134.0 

6.0 

NGC  7026 . 

[N  ii] 

51 

52 

3.00 

0.17 

0.67 

6.00 

99.7 

6.0 

27.5 

[O  m] 

53 

54 

3.00 

0.17 

0.67 

6.00 

99.7 

6.0 

24.0 

NGC  7354 . 

[O  m] 

55 

56 

0.50 

0.40 

0.66 

0.75 

91.7 

32.0 

23.0 

Ho 

ZK 

ZK 

1.00 

0.20 

0.94 

20.0 

Notes. — The  columns  represent  PN  name,  species,  the  number  of  the  figure  containing  the  data  and  the  model  fit,  respectively  ( ZK  refers  to  the  fit  obtained  by  Zhang  & 
Kwok  1 998),  a  ( latitudinal  density),  (3  ( polar  to  equatorial  density  ratio),  A  r  (axial  ratio),  Iq  (ionization  parameter),  i  (inclination  angle  of  major  axis),  O'  ( position  angle  of 
major  axis  on  sky),  and  vq  (equatorial  expansion  velocity).  In  all  cases,  the  parameter  7  (exponent  of  radial  density  power  law  within  the  rim)  is  set  to  0. 


lengths  of  the  nebular  major  and  minor  axes.  Because  the  ob¬ 
served  nebular  geometry  is  a  strong  function  of  inclination  angle, 
the  incorporation  of  velocity  data  into  a  model  significantly  con¬ 
strains  several  parameters,  including  the  axial  ratio.  Of  the  nine 
PNe  that  we  have  in  common  with  the  analysis  performed  by 
Zhang  &  Kwok  (1998),  seven  had  major  axis  lengths  that  dis¬ 
agree  with  our  results  by  more  than  the  1  er  uncertainty  ( 1 0%).  Of 
the  eight  nonspherical  PNe  (1C  3568  is  spherical,  and  therefore 
its  inclination  value  is  irrelevant)  seven  have  inclination  angles 
that  exhibit  disagreement  by  more  than  1  a  (20°).  And  only  four 
of  these  eight  PNe  have  axial  ratios  that  are  in  agreement  within 
1  a  (0.15).  These  results  illustrate  the  importance  of  using  the 
velocity  data  to  determine  accurate  nebular  parameters  and  will 
have  important  implications  in  deriving  direct  expansion  distances 
to  these  objects. 

Although  the  quality  of  our  model  results  is  generally  a  good 
match  to  the  data,  some  of  the  kinematic  data  show  line  splittings 
that  the  model  cannot  reproduce.  For  example,  consider  Fig¬ 


ures  3  and  4.  The  data  for  the  “A”  slit  position  in  Figure  3  clearly 
shows  the  [O  in]  emission  at  two  discrete  velocities,  as  if  a  re¬ 
ceding  and  an  approaching  wall  of  gas  were  present.  However, 
the  kinematic  model  along  this  slit  position  in  Figure  4  shows  the 
line  being  split  in  the  spatial  dimension  (i.e.,  the  emission  peaks 
occur  at  discrete  positions).  We  were  unable  to  remove  this  effect 
despite  considering  a  variety  of  model  parameters,  and  we  con¬ 
clude  that  the  discrepancy  is  due  to  a  deficiency  in  the  model, 
perhaps  in  the  way  that  the  gas  distribution  within  the  thin  neb¬ 
ular  shell  is  treated.  The  discrepancy  between  the  model  and  data 
also  appears  to  be  present  in  some  of  the  [O  m]  images/spectra 
but  surprisingly  is  not  common  in  [N  n]  data  sets.  This  may  be  a 
clue  to  eventually  tracking  down  the  source  of  the  disagreement. 

4.  SUMMARY 

In  this  compilation,  we  present  a  catalog  of  images  from  the 
WFPC2  on  board  HST,  and  long-slit  spectra  at  various  positions 
across  a  sample  of  PNe  in  the  light  of  [N  n]  and  [O  m].  Many/most 
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of  the  PNe  in  this  compilation  are  well  represented  by  ellipsoidal 
geometries.  For  these  targets,  we  fit  models  based  on  the  PES 
model  by  Aaquist  &  Kwok  (1996)  and  Zhang  &  Kwok  (1998), 
and  modified  to  predict  the  nebular  kinematics,  as  well  as  spatial 
morphologies.  A  comparison  between  our  best-fit  parameters  and 
the  above-mentioned  studies  shows  differences,  mainly  resulting 
from  inclusion  of  the  kinematic  modeling  and  resulting  parameter 
constraints  that  are  not  provided  by  the  image  data  alone.  The 
inclination  angle  and  axial  ratio  of  the  nebula  are  strongly  af¬ 
fected  by  inclusion  of  kinematic  data,  and  most  of  our  estimates 
are  in  complete  disagreement  with  estimates  based  only  on  im¬ 
age  data.  The  inclination  angle  of  the  nebular  major  axis  relative 
to  the  line  of  sight  is  of  special  interest  to  us,  since  our  ultimate 
goal  is  to  infer  expansion  distances  to  these  PNe.  In  addition,  the 
substantially  increased  spatial  resolution  offered  by  the  Hubble 
Space  Telescope  as  compared  to  ground-based  observations  al¬ 
lows  us  to  completely  resolve  the  nebular  shells.  Not  surpris¬ 
ingly,  this  has  a  major  impact  on  estimates  of  the  nebular  shell 
thickness.  The  precision  of  the  HST  images  and  the  addition 
of  kinematic  data  result  in  better  estimates  of  nebular  paramet¬ 
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ers,  which  will  enable  us  to  determine  accurate  distances  to  the 
selected  PNe. 
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